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Abstract 



Collision cross sections and rate constants between two ground- state oxygen 
molecules are investigated theoretically at translational energies below ~ IK 
and in zero magnetic field. We present calculations for elastic and spin- chang- 
ing inelastic collision rates for different isotopic combinations of oxygen atoms 
as a prelude to understanding their collisional stability in ultracold magnetic 
traps. A numerical analysis has been made in the framework of a rigid- rotor 
model that accounts fully for the singlet, triplet, and quintet potential energy 
surfaces in this system. The results offer insights into the effectiveness of evap- 
orative cooling and the properties of molecular Bose- Einstein condensates, as 
well as estimates of collisional lifetimes in magnetic traps. Specifically, 17 O2 
looks like a good candidate for ultracold studies, while 16 O2 is unlikely to 
survive evaporative cooling. Since O2 is representative of a wide class of 
molecules that are paramagnetic in their ground state we conclude that many 
molecules can be successfully magnetically trapped at ultralow temperatures. 
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I. INTRODUCTION 



A. Background 

Following the enormous successes of lowering the temperature of atoms to the sub-mK 
regime, experimental attention is now turning to producing ultracold molecular samples. 
The "first generation" of cold molecule experiments has now demonstrated the efficient 
production of cold samples by a variety of techniques, including photoassociation of ultracold 
atoms |J, counter-rotating supersonic jets 0, Stark slowing 0, and buffer-gas cooling [|J. 
The latter two have yielded trapped samples that are cold in rotational, vibrational and 
translational degrees of freedom, although translational temperatures are still in the 0.1-1 
K range. 

The next generation of experiments will seek colder, denser samples. One novel option for 
cooling molecules further is the optical cooling strategy described in Ref . || . Alternatively, 
we consider in this paper the evaporative cooling of paramagnetic molecules in a static 
magnetic trap, following the strategies that have been employed to produce ultracold atoms. 
The central issue to the success of this method is that the rate of elastic, rethermalizing 
collisions far exceeds the rate of lossy inelastic collisions that produce untrapped, strong- 
field-seeking states. In terms of scattering rate coefficients, this criterion is usually written 
K e \ > 10 2 A"i oss . A large ratio of K e \/Ki oss is also vital for the stability of the trapped gas 
once it is cold. A main objective of this paper is to demonstrate that molecules with nonzero 
spin in their lowest energy state will be quite stable at ultralow temperatures. 

We may reasonably assert that the elastic rate constants for neutral molecules have 
roughly the same magnitudes as those for neutral atoms, K e \ ~ 10~ 12 — 10~ 10 cm 3 /sec at 
low energies, barring unfortunately placed zeros in their s-wave scattering cross sections. 
Indeed, our calculations yield elastic rates of this magnitude. The spin- state- changing 
rate constants are, however, completely unknown for molecules at ultralow temperatures. 
To rectify this situation the present paper presents pilot calculations for cold collisions of 
molecular oxygen. This work is a logical next step following Refs. |6|,[7]], which considered 
the interaction of molecular oxygen a the helium buffer gas. 

Spin- state- changing cold collisions of O2 molecules are driven by at least three distinct 
physical processes, two of which are already familiar from ultracold atom physics: i) Spin- 
exchange collisions, which typically lead to unacceptably large loss rates for both atoms and 
molecules; ii) Spin-spin magnetic dipolar interactions, which are typically small in either 
case; and iii) Spin-rotation interactions, unique to molecules, wherein electronic spins are 
influenced by their coupling to rotational motion, which is in turn dependent on torques 
exerted by the anisotropic potential energy surface(PES) between the molecules. 

Spin-exchange collisions can be avoided in cold molecule collisions, as in cold atom col- 
lisions, by preparing the molecules in their "stretched" states, with mechanical rotation, 
electronic spin, and nuclear spin (if any) all aligned along a common laboratory-fixed axis. 
We will therefore confine ourselves to this circumstance. By far the leading contribution 
to the rate constant for state- changing, lossy collisions (Ki oss ) is then the spin-rotation 
coupling, as shown below in detailed calculations. Indeed, when an exothermic exit channel 
is available, this coupling can yield loss rates comparable in magnitude to spin-exchange 
rates, i.e., comparable to elastic collision rates. This is the case for the 16 C>2 molecule. The 



2 



spin- changing rate is, however, strongly suppressed when the only allowed exit channels are 
degenerate in energy with the incident channel and when the collision energy lies below a 
characteristic energy Eq. In the stretched state of 17 C>2 this is indeed the case, since chang- 
ing the molecular spin at low energy requires boosting the partial- wave angular momentum 
from I = to I = 2. These collisions are therefore suppressed by the Wigner threshold law 
when E < Eq, where Eq is the height of the I = 2 centrifugal barrier. For 17 02 the barrier is 
roughly Eq ~ 0.013 K, not far below the temperature that buffer-gas cooling can take these 
molecules. 

The potential for disaster in molecule cold collisions is far greater than in atom cold col- 
lisions. For example, hyperfme interactions are more complex, and include rotation-nuclear 
spin couplings that can invert spins. These are, however, expected to have minimal impact 
on the stretched-state molecules. Also of potential significance are spin-vibration couplings, 
which we disregard in O2 owing to the extremely large vibrational excitation energy of 
ground state O2 molecules compared to the energy available to excite them. The vibrational 
degrees of freedom remain to be fully explored at ultralow temperatures, but it is expected 
that vibrational quenching (also an exothermic process) can occur with appreciable rates ||. 
Finally, polar molecules are susceptible to particularly strong long-range anisotropies. While 
this is not of direct relevance to molecular oxygen, it can be devastating to the electrostatic 
trapping of polar molecules || , and potentially dangerous for magnetically trapped molecules 
when the electric and molecular dipole moments are coupled. 



B. Oxygen molecules- General considerations 

The importance of molecular oxygen as a potential candidate for cooling and trapping 
experiments has been emphasized elsewhere PJT0|j. We will here consider O2 molecules 
that have been cooled to temperatures below IK. We will furthermore assume that these 
molecules have relaxed to their electronic 3 X~ ground state, and v — vibrational ground 
state. What remain are the rotational and spin degrees of freedom that influence evaporative 
cooling. In a typical magnetostatic trap the molecules can be confined provided that they 
are in a weak-field-seeking state, i.e., a state whose energy rises with magnetic field. 

The Zeeman diagram of O2 is reproduced in Figure 1. Nuclear exchange symmetry 
declares that homonuclear oxygen isotopomers can have only even or only odd values of the 
nuclear rotation quantum number AT [12|. For 16 C>2 (or 18 C>2), which has identical spin- zero 
nuclei, only odd values of N are allowed. For the isotopomer 17 C>2, where each nucleus has 
spin I N = 5/2, the allowed iV levels depend on the total nuclear spin I = I N a + Inb- The 
molecules must have even values of N for odd values of J, and vice-versa. Thus the lowest 
energy weak-field seeking states are \N,JMj) = 1 1, 22) and 1 1 , 21) for 16 C>2 or 18 02, and 
\N, JMj) = |0, 11) for 17 C>2. These states are indicated by heavy lines in the figure. 

Figure 1 illustrates the essential difference between the even-iV and odd-iV manifolds in 
O2, from the standpoint of inelastic collisions. Namely, the trapped states with J = 2 in the 
N — 1 manifold can decay exothermically to the untrapped J = states. By contrast, the 
N = trapped state at low energy and low magnetic field can only change its spin projection 
Mj to other states that are nearly degenerate in energy. This difference proves crucial in 
strongly suppressing spin-rotation collisions in 17 C>2 relative to 16 02- This has already been 
discussed for cold collisions of O2 with helium atoms PJ7I1; the situation is similar when 
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the molecules collide with each other. As in Refs. we carry out calculations in zero 

magnetic field. 

In this paper we ignore the role of nuclear spin, hence of hyperfine structure, in the 17 C>2 
isotopomer. This is justified by considering the molecules to be in their stretched states of 
N = 0, with Mj = J = 1 and Mi = 1 = 5. The nuclear spin degrees of freedom are then 
frozen out. In principle the weak nuclear spin-rotation coupling would influence the nuclear 
spins, resulting inF = 6^F = 5,4 transitions, but these couplings are an order- of- 
magnitude smaller than the spin- rotation couplings we are already considering. Moreover, 
measurements of the 17 O2 microwave spectrum ]TT| reveal that its hyperfine structure is 
inverted in its N = state, requiring that 3.7mK of kinetic energy be supplied to change F. 
Thus at ultralow temperatures hyperfine- state- changing collisions are rigorously forbidden. 



II. MODEL 
A. Hamiltonian 

Our investigation of molecular collisions follows the model of diatom-diatom scattering 
originally due to Hl3|,[l4|, but modified to incorporate the electronic spin of oxygen molecules. 



The 02( 3 E ff ) — 02{ 3 ^ g ) dimer has a spin-dependent intermolecular potential, namely three 
potential surfaces exist corresponding to singlet(5' = 0), triplet^ = 1), and quintet (S = 2) 



states of total electronic spin S |15[| . The complete Hamiltonian for the collision process can 
be written 

H = T A + T B + V s + V lr + V dd + H rfs , (1) 

where Tj is the translational kinetic energy of molecule i; V s is the short-range exchange 
interaction; Vi r is the long-range potential consisting of dispersion and electric quadrupole- 
quadrupole interactions; Vdd is the electronic spin-spin dipolar interaction; and H r f s is the 
Hamiltonian for the rotational fine structure of the two separate oxygen molecules. The 
hyperfine interaction will be disregarded for now, as was discussed above. The short-range 
potential can be written as a mean interaction plus exchange corrections, following Ref. pi : 



V s (R,ua,u}b, Sa, S b ) = V av {R,u A ,u B ) - 2V r ex (R, u A , u B )S A ■ S B , (2) 

where to = (9, 0) are the polar angles of molecules A and B respectively, R(i?, 0) describes 
the radius vector between the center-of-mass of the molecules in the lab- fixed coordinate 
frame and 



E fL A ,L B AR)(-i) LB ~ L y(2L A + 1)(2L B + 1)K L (^,^)-C L (6). (3) 

La,Lb,L 

Here (u A , uj b ) = [Cm a a {uj a ) ® C^ b {ijJb)}m an d Cm are reduced spherical harmonics. It 
should be said that the expansion (||) is identical to that in [16] but written in terms of re- 
duced spherical harmonics. The spin-dependent Heisenberg exchange term V ex is expanded 
similarly, but with a different expansion coefficient gL A ,L B ,h{R)- The expressions for ex- 



pansion coefficients fL A ,L B ,L(R) and gL A ,L B ,L(R) were obtained in the work of |15| , and the 
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quadrupole- quadrupole interaction has a similar form. The G% dispersion coefficients were 



calculated in |17| in the body- fixed frame. To unify our treatment we recast the anisotropic 



Cq coefficients in terms of the same angular basis as the exchange potential: 



Pl a ,L b ,L ( 1 \L 

La,Lb,L 



E ^^(-l) LB - LA ^/(2L A + l)(2LB + l)K L (u A ,u B yC^Q), (4) 



where 



PLa ' Lb ' L ~ V (2L A + 1)(2L B + 1) J£ Mb d ^M A ,Un [ m A . aL n ) ( 5 ) 




The connection between our coefficients dB A ,L B MAMB and coefficients from [JlT] is in Ap- 
pendix A. Figure 2 shows a slice through the potential energy surface for the singlet states 
of the [02( 3 Sg )]2 dimer for the "H" -geometry of the two molecules, which passes through 
the global minimum of the PES. Also shown are contributions to this potential from various 
components with different (L, L A , Lb)- We can see that this potential has a very strong 
anisotropy. Notice that the isotropic contribution, with L A = Lb = L = accounts for less 
than half of the total well depth. 

As in the case of ultracold atoms, details of ultracold molecular collisions depend ex- 
tremely sensitively on details of PES. Thus eventually the PES must be fine- tuned using 
key experimental data to provide complete quantitative results |TB,n5|. Nevertheless, the 



order-of- magnitude of the rate constants, and their general trends, already emerge clearly 
in the present model. 



The intermolecular spin-spin (magnetic dipole) interaction has the form ]16 



V M (R, S A , S B ) = -^^C 2 (6)-[S A ® S B ] 2 , (6) 

where R = (R, G), g e = 2.0023 and \xb is the Bohr magneton. For alkali atoms this is known 
to be a fairly weak contribution to spin-changing collisions, a conclusion that we find holds 
for molecules as well. Finally, in our model we take into account the molecular fine structure, 
which arises from the molecular rotation and spin-rotation coupling and is diagonal in our 
total-spin basis at large R. For 16 02 we use the fine- structure constants determined in [|20 
and for 17 O2 we employ those determined in |21 



We express the Hamiltonian in a basis of total angular momentum, 

J = J + l 

J = Ji + J 2 (7) 
Ji = Si + N i; 

in the terms of each molecule's mechanical rotation (Ni), its electronic spin (Si), its total spin 
(Ji), the combined spin for two molecules together (J), and the partial wave representing 
the rotation of the molecule about the center of mass (1). In zero magnetic field both J and 
its laboratory-fixed projection M. are rigorously conserved. 
In this basis we can present our wave function as: 
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U l,J,Ji,J 2 ,Ni,N 2 ,Si,S 2 

Ij,M\l,J,Jt,J2,Ni,N 2 ,Si,S-2{®i u Ai U B, &A, CTb), (8) 

where oa,b represent the electronic spin coordinates for molecule A,B. The coupled angular 
momentum basis functions are defined by suitable tensor products: 

Ij,M;i,J,Ji,J2,Ni,N 2 ,Si,s 2 (@,UA,UB,crA,VB) = [Pmj(ua,ub,(7a,(7b) <§> (6)]^, (9) 
P Mj (uJ A ,uj B} a A ,a B ) = \Tm 3i (uj a ,ct a )® (uj B ,a B )} J Mj} (10) 

T£ Jf (u,<r) = [Y$ Ni (w) ® , (11) 

where Y is a spherical harmonic and \ is a spinor wave function. 

Because target and projectile are identical bosons, we must take into account the sym- 
metry of the wave function under exchange. To this end we construct symmetrized angular 
momentum functions from (|9| ): 

T s _ ^l,J,Ji,J2,N 1 ,N 2 ,S 1 ,S2 + {-^) J+Jl+J2+l Il,J,J 2 ,Ji,N 2 ,N 1 ,S 2 ,S 1 / 10 x 
1 l,J,J 1 ,J 2 ,N 1 ,N 2 ,S u S 2 - r— \ LZ ) 

y 2(1 + dj^djf^dsiSi) 

We have omitted the conserved quantum numbers J , M. in this expression. 

To calculate the reduced matrix elements in our basis we recouple the angular part of 
the potential in terms of tensor elements: 

L(,.,. ,., \rc. . g_17 _ . I\ TJ\ t3(^\l-l a -j b J L A L B 7 



K^ A ,u B )[S A • SbF = E \J[Ja][Jb](-1) L - La - Jb [ j b j a l] ( 13 ) 
in terms of the spherical tensors T(f (uj a ,<j a ) = [Y^f ® [Xmo ] 7 ]m 7 , with 7 = for V av 

J A ''A -'A 

and 7 = 1 for the exchange part of the interaction (|2|), and [Q] = (2Q + 1). 

Using our expansion of the inter- molecular potential (^ and (^), the wave function @ 
and taking into account the Wigner - Eckart theorem , 

< JM\V 3 + V disp + V dd \J'M' >= 5j J/ 5 MM/ < J\\V S + V disp + V dd \\j' >, (14) 

we can present the reduced matrix elements for the V s + Vi r part as: 

< {[J 1 (N 1 S 1 )J 2 (N 2 S 2 )]Jl}J\\(K L ■ C L )[S A • S B r\\{[J[(N[S[)4(N^S' 2 )]J'l'}J > 

= (-1)*+'-^+*+* ( [i'} [f] [i] [j] [j> 2 \ [j,] [ j 2 ] [n[] m m [n 2 ] 

x / / V L \ ( I L l' \ ( N x L A N[ \ ( N 2 L B N 2 \ 

lj'jj'jloooMo M o j 



(15) 



C t t 1 f J J' L j r J\ J[ J A j r J 2 J 2 J B 

E (-1) Jb [Ja][Jb] { j j J Ji J[ Ja M N[ L a \{ N 2 N> L b 

A } { J2 J' 2 Jb ) { S 1 S[ 7 J { S 2 S' 2 7 J 

S A (S A + 1)(2S A + 1)S B (S B + 1)(2S B + VjYds^Ss^, 



x 

JaJb 



6 



and for the Vm part as: 



< {[Ji(N 1 S 1 )J 2 (N 2 S 2 )]Jl}J\\(C 2 • [S A ® S B ] 2 )||{[Ji(NiS' 1 )J , 2 (N , 2 S^)]J'l , } t 7 > 

= (-i) j+j,+jVi+iV2+Ji+j2 ([n[Ji[/][j] [fMuw^ 2 [y'jj] ( 16 ) 

J J' 2 

Jl J 2 1 

J{ 1 

x y/S A (S A + 1)(25 A + l)S B {S B + 1)(2S'b + 1)^^^^, 

The reduced matrix elements of our potential between the states defined by (|j) and 
using (O) are 




< nhW'Wrfj'j' >= - ^[WjJ > +(-i) J+Ji+j2+; < ^HgbMj > 

^(1 + Sj lJa 6 Nl N 2 )(l + fyj^iv^ 

x +l } 2 , (17) 

where r\ stands for the rest of the quantum numbers from our wave function (|S]). The coupling 
matrix element therefore vanishes between channels with different parity (— l^i+^+z . 

Figure 3 shows a partial set of adiabatic potential curves for 16 2 in the case of J = 0. 
To generate this figure we include rotational channels N = 1, 3, 5 and even partial waves 
I = — 6, which already imply 100 channels in this case. The strong anisotropy in the PES 
is here manifested mainly in a set of strongly avoided crossings near R = 8 a.u. For R > 8 
a.u. the PES strongly favors a collinear configuration of the pair of molecules, while for 
R < 8 a.u. it strongly favors a parallel, "H"-shaped configuration. This figure stresses the 
importance of higher-lying rotational states in determining the details of scattering even at 
ultracold energies. However, for 16 2 we have chosen to compute cross sections just for the 
case I = — 10, N = 1 because these calculations already reveal very large spin-changing 
rates. Higher-lying channels will influence the details, but are unlikely to suppress losses. 
The total number of channels, including all values of J , is then 212. For 17 2 , by contrast, 
we have computed cross sections for I = — 10 and iV = 0, 2, to verify that higher-lying 
rotational states do not upset the observed suppression of loss rates. In this case the total 
number of channels considered is therefore 836. 



B. Evaluating cross sections 



We solve the coupled-channel equations using a log-derivative propagator method |22 
to determine scattering matrices. Since we assume zero magnetic field the total angular 
momentum J is a good quantum number and the results are independent of the laboratory 
projection Ai of total angular momentum. 

For magnetic trapping the molecular quantum numbers of interest are naturally the 
magnetic quantum numbers. Therefore we need to know the state- to- state cross sections 
in the \NxN 2 J\ J 2 , M Jx Mj^ > basis. The scattering matrices are readily converted to this 
basis: 
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< N^JihMjiMjJMilSlNlN^JzM^M'j/Ml >= 

E < JiMj x J 2 Mj 2 \JMj >< J'M'j^M'j^M'^ > 
jj' 

x E < JMjlMi\JMj >< JMj\J'M'jl'M{ > (18) 
J 

x < {[J 1 (iV 1 ^ 1 )J 2 (iV 2i 5 2 )]Ji} t 7^|<S( l 7)|{[J((iV^04(iV^)]JY} l 7^ > 

For notational simplicity we define the index a = (NiN 2 JiJ 2 ) in the following. Then the 
complete symmetrized wave function in the limit of large R is given in [ I3| : 



(exp(ik a ■ R)T^l h (ua, oa)Tm 3% (wb, cr B ) + exp(-ik a ■ R)T^ (u B , (Tb)^ {uj a , <?a) 

exp(tk' a R) fj[Mi J>Mi {R) + fj>M> JjM* (-R) j, j, 

+ E p 1 7k 1 t m,, ( u a, va)T^, (u b , ct b ) (19) 

a'M'M' K V2 J i J 2 

J l J 2 

where fji m> ji m' is the channel-dependent scattering amplitude. 

1 Ji 2 j 2 

Using the definition of our wave function (B) and transforming it into the | J\ J 2 Mj x Mj 2 > 



basis we can get the asymptotic form of the wave function in terms of the S- matrix [p^| : 



E < JiMj, J 2 M j2 1 JMj >< JMjlM t \JMj > y/1 + Sj iJ2 6 NiN2 5 Mji mj 2 

J,J,l,Mi 

x^^i'Y^ik) E 4= E (Sj^jrSjj^expi-iik^R-l'Tr^) (20) 
-exp(i(k' a - R-l'ir/2) < {[J X J 2 \ Jl}JM\S{J)\{[J[ J' 2 \ J'l'}JM >) 

xIjM-Vj'^j'Jy®, UA, UB, <7b), 



By comparing (|T^) with (PD|), we obtain the expression for the scattering amplitude: 

Air 



fj[M' h J' 2 M' j2 {R) + fj^M'j^M'j^-R) 



k a k' a 



X E \A + Sji^N^M^M^J^ + h[J^N[N^M'j M'j ^,(4)^ (21) 

x < N^JxMj^MjJM^S-I^N'zJiM^MjjMl > Y&,(R), 
where in the symmetrized separate- molecule basis the scattering matrix is given by 

< NiNzJiJzM^MjJMilSlN^J^M'^M'j/Ml >= 

\ 1 + 5j, J 7 5]\f, M> - . 

v " ' E < ■hMj 1 J 2 Mj 2 \JMj >< J'M'j\J' x M'j x J' 2 M'j 2 > 



\J 1 + S J x J 2 5 Nl N2 S Mjl Mj 2 



JJ' 



1 + Sj' J'Sn' N' 

, r r E < JMjlM^JMj >< JMj\J'M'jl'M[ > (22) 

x < {[J 1 (iV 1 5 1 )J 2 (iV 2 5 2 )]J/} l 7A^|5( l 7)|{[J((iV^0J 2 (iV 2 5 2 )]J7 / } l 7^ > . 



Symmetrization in the | J\JiMj x Mj^) basis requires that J\ > J2 and that M\ > M2 when 

To obtain the scattering cross section we must integrate over the angular coordinates of 
the scattered wave. But, for undistinguishable final spin states we restrict the integral over 
half space (/ dQ = 2n) to avoid double counting p3| . 

The total state- to- state cross section of interest for spin-rotational excitation and re- 
laxation phenomena can be obtained using ([H]) from the S matrix: 

@(NiN 2 )JiJ 2 Mj 1 Mj 2 ^(N' 1 N 2 )J[J 2 M' J M'j = 72 

K N 1 N 2 .hJ2 

x J2 I < {NxN 2 )J l MjJ 2 Mj 2 lM l \S-l\{N[N' 2 )J[M' Ji J , 2 M'jM' l > | 2 , (23) 

lMil'M[ 

where 

kN 1 N 2 JiJ 2 — ~~ E NlJl - E N , 2 j 2 )) 1/2 (24) 

is the channel wavenumber and E^ 12 j 12 are fine structure energy levels. In this expression 
we assume an average over all incident directions, as in ||. Finally, state- to- state rate 
coefficients are given by 

K(N 1 N 2 )J 1 J 2 M Jl Mj 2 ^(N[N 2 )J[J' 2 M' Ji M' j2 = V( Nl N 2 ) J x J 2 ^(N X N 2 ) J x J 2 M Jl Mj^N^J^M'^ A/^ , (25) 

where V[n 1 n 2 )j 1 j 2 is the relative velocity of the collision partners before the collision. 

III. RESULTS 

This paper considers the scattering problem for the homonuclear species 16 O2 and 17 02- 
Ref. H[7[ speculated that buffer- gas cooling of O2 by Helium should be possible, thus 
lowering the molecules to typical temperatures ~ 0.3K To further cool the gas by evaporative 
cooling requires favorabale collision rates for collision energies E < IK. We will limit our 
detailed calculations to this case. We will see that the cooling of 11 02 could be quite efficient, 
while it is probably impossible for 16 2 - 

A. 17 C>2- elastic scattering. 

Since 17 O2 is the most promising candidate for evaporative cooling, we devote our atten- 
tion to this isotopomer. We focus on the \N1N2J\J2, Mj x Mj 2 >= 1 001 1, 11 > state which is 
the lowest- lying trappable state for the even N- manifold (Fig.l). 

For identical bosons only even I- partial waves contribute to the cross sections for the 
|0011, 11 > state. A first important point of our calculations is to determine the number of 
partial waves that contribute to the cross section in the energy region up to IK and how 
many molecular rotational levels should be taken into account. In principle many partial 
waves are coupled together by the very anisotropic potential, but higher partial waves are 
suppressed at a low energy. Figure 4 illustrates the elastic cross sections for different valves 
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of the highest partial wave included. This Figure shows that it is enough to include just 
1=0-4 partial waves for the qualitative description of the cross section in the region up to 
~ 0.2K, and that the partial waves 1=0-10 are sufficient in the region up to ~ IK. For 
all our calculations for 17 2 we considered just the two lowest rotational levels N = 0,2. 
Including only the N = rotational level allows the molecules to explore only the isotropic 
part of their PES. Thus N = 2 states must be included at least. The influence of higher 
rotational levels are found to be small in test calculations, although they impact details of 
the resonance structure. The calculations thus include 836 channels. 

Particularly striking in Figure 4 is the strong difference in the cross sections when l max = 2 
as opposed to l max = 4, 6, 8, 10. This is caused by the strong anisotropy of the potential fl5( 
and mathematically this means that for the case l max = 2 we take into account only a small 
number of expansion functions (L, La, Lb) in (|3|). Although, these few functions represent 
"most" of the potential, even small changes in potential can change the behavior of the 
cross section dramatically near zero energy TP| . Figure 5 shows the elastic cross sections 
for 1 00 11, 11 > collisions and the contribution from different partial waves. We can see that 
the partial waves I = 8 and I = 10 contribute significantly only above ~ 0.5K. 

Figure 5 also exhibits dozens of resonances for molecular- molecular collisions below IK, 
arising from the enormous number of internal molecular states JF|. Although we do not 
assign quantum numbers to the resonant states here, we expect them to be of two basic 
types, as discussed in [J/J: i) coupled- channel shape resonances; ii) "rotational Feshbach" 
resonances that change the value of N of one or more molecules. This last type of resonance 
can be extremely long- lived, owing to the difficulty of both molecules returning to their 
rotationless state in a collision. Indeed, one such resonance at ~ 0.013-fT has a width 50fiK. 
We will return to this subject in a future publication. 

Figure 5 shows that the elastic cross section has a very large value near zero energy, 
corresponding to a scattering length a = 270a. u. in the present model. It is therefore 
possible that there may be an s- wave bound state in the region of negative energies near 
the threshold of the channel. In this case the cross section should have ~ 1/(E + \e\) 
dependence [24| on energy E and on the energy of the bound state e. 

We have included the dipole- dipole interaction in these calculations. However the role 
of this interaction is very small in general, influencing the cross section at the 1% level and 
shifting resonance positions slightly. Likewise, this interaction is only a small perturbation 
to inelastic scattering. 



B. O2- prospects for evaporative cooling. 

For the | N1N2J1J2, Mj x Mj 2 >= 1 001 1,11 > state of interest to trapping experiments, 
Ni and Jj are conserved at low energy, since the next energetically available state (with 
iV = 2) is 11.18K higher in energy (Fig.l). Thus the only possible final states are those 
which differ from the initial one in their projections Mj 1 and Mj 2 . To accomplish such 
a transition therefore requires that the angular momentum be carried away in the orbital 
angular momentum I. Furthermore, the collisions that originate in s- wave channels will be 
suppressed at energies below the centrifugal barrier of the d- wave exit channel. Using an 
effective Cg coefficient, Cg from |I7| this energy for a partial wave I can be approximated 
as 
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6Cq m 
H 2 l(l + 1) 



(26) 



For 2 , Cl n = 80.5 a.u. and the d-wave threshold energy is 0.013/T. 

The main aim of this paper, as previously discussed, is to compare the elastic and loss 
rate constants. Figure 6a shows these rates calculated according to Eqn.(pBl). Away from 
resonances, in the energy range up to E w 0. 013.fr the loss rate constant is indeed strongly 
suppressed. A detailed examination of the final states contributing to this loss reveals general 
features which are the similar to those for He — 17 O2 scattering. Namely, elastic scattering, 
which does not change either Mj 1 or Mj 2 , is the most probable result of a collision. The 
next most likely processes are those for which the final Mj differs from the initial one by 1 
or 2 and the rates for these processes are smaller than the elastic rate by one -two orders 
of magnitude. The processes for which the final Mj differs from the initial one by 3 or 4 
have rates smaller than for elastic scattering by 2-4 order of magnitude. For energies above 
« 0.0137^ inelastic processes become more probable, with rates only about 7- 10 times 
smaller than elastic rates. 

The thermally averaged elastic and loss rates are relevant to the experimental situation. 
If we assume the velocity distribution is Maxwellian characterized by a kinetic temperature 
T we can calculate the thermally- averaged rate constant as: 

*< r H^F) 1/2 (^f iW*^ (27) 

To do this averaging we extrapolate the cross sections to energies > IK using their values 
at E = IK. 

Figure 6b shows these thermally averaged elastic and loss rates. For the cooling to 
be efficient the rate of elastic collisions K e \ must exceed the rate of spin- changing, lossy 
collisions K\ oss by at least two orders of magnitude p^ . For 17 2 , below T « 0.01-fT this 
condition is fullfilled. However, there is a "relatively dangerous temperature range" above 
« O.Oli^ where K e i/Ki oss m 7—10. By comparison, consider the equivalent ratio for He — 02 
scattering, as discussed in Ref . fjj . The fact that K e i / Ki oss is not so large as for He — O2 
collisions originates from the stronger anisotropy and the deeper PES for the O2 — O2 system. 
It remains to be seen if the loss rates are sufficiently low to evaporatively cool from buffer- 
gas temperatures, ss 0.3K, down to T < O.Oli^, where cooling should be quite efficient. 

C. 16 o 2 

It is a different situation for 16 02 molecules from the point of view of comparing elastic 
and inelastic cross sections. The general behavior of the elastic cross section for different 
channels is similar to that of 17 O2 and has the same order of magnitude except in the energy 
region near zero which is very sensitive to the details of the potential and the reduced mass. 
With the present PES the 16 02 scattering length is 28a.u. 

Figure 7 shows the elastic and all the inelastic cross sections for the trapped state 
|iViJV 2 J r i^2,Afj 1 Mj a >= |1122, 22 >. The total number of inelastic channels is 25. When 
the final states are 1 1120, 20 > , 1 1100, 00 >, 1 1110, 10 > and so on, i.e when at least one 
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of the molecules changes to the J = state, the collision is superelastic. It is well known 
that for a superelastic channel there is "cr ~ 1/v" threshold law. Thus at a low energy 
there is a substantial loss of molecules from the 1 1 122, 22 > state. The same result holds for 
1 1122, 11 > state which is also, of course, susceptible to spin- exchange. Thus 16 2 is clearly 
unstable against collisional losses in a magnetic trap, in sharp contrast to 17 O2. The cross 
sections for the stretched states of 16 02 molecules that we are interested in have a smooth 
structure in the energy region up to IK, implying either a lack of resonances in this region or 
their large widths. From Figure 7 we can see just one sharp Feshbach- type resonance near 
0.87^ belonging to an I — 10 bound state. Although there may be other weak resonances, 
it was not our aim to find all the resonances and identify their nature in this article. 



IV. CONCLUSION 

In this article we theoretically investigated ground state diatom- diatom collisions in the 
energy range up to IK taking different isotopomers of oxygen molecules as a prototype. 
The main point of our investigation was to estimate the ratio of elastic and inelastic rate 
constants. The influence of the rotational degrees of freedom is crucial in determining 
this ratio. In the case of the odd N- manifold, it is probably impossible to satisfy the 
criterion K e \ > 10 2 K loss for stretched states in any energy region, because of both the strong 
anisotropy of the PES and the existence of the superelastic channels. In the case of the even 
N- manifold, namely N = 0, this criterion can be fullfilled because for the stretched state 
1 00 11, 11 > there are no superelastic channels. 

Even though the required ratio of K e \/K\ oss is not quite met at buffer gas temperatures, 
it is worth remembering that the buffer-gas procedure typically produces a far larger sample 
of trapped molecules than do the laser cooling experiments on which evaporative cooling 
is usually applied. Thus it is possible that a larger loss rate could be sustained without 
harming the overall yield of molecules at ultralow temperatures. Detailed rate-equation 
simulations of the cooling process are therefore required, an item to which we will turn our 
attention in the future. 

Equally importantly, once the molecules have in fact been cooled to /xK temperatures, 
out results imply that the lossy collision rates have diminished into insignificance, falling to 
levels well below 10~ 14 cm 3 /sec. This in turn implies that ultracold spin-polarized 17 O2 gases, 
like their atomic counterparts, are experimentally stable and should allow the production of 
novel Bose-Einstein condensates. 

Beyond the immediate results for our particular model of 17 02, the present results have 
broad implications for many paramagnetic molecular species. Namely, the characteristic 
suppression of loss rates below the d-wave centrifugal barrier should be a generic feature for 
molecules where no superelastic fine-structure-changing processes exist. It is also important 
to assess in detail the influence of the anisotropy of the PES. For this purpose further 
investigations are necessary. 

This work was supported by the National Science Foundation. 

Appendix: The definition of the coefficients. 



Using the dispersion contribution part from JT7| and expanding this part in the terms of 



spherical reduced harmonics (m in the lab-fixed coordinate frame we defined the connection 



12 



between our coefficients cIl a ,l b ,MaMb an d coefficients A,B,C from fl7 



4,o,o,o = + 8B + 8C); 4,o,2,o = 2A-AB + 2C 

4,o,o,o = 4,0,2,0 = ^(2A + 2B- AC) 
4,-i, 2 ,i = 4,i,2,-i = \{^A - 8B + AC) 

4,-2,2,2 = 4,2,2,-2 = \{2A-AB + 2C) 
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FIGURES 



FIG. 1. The lowest-energy Zeeman levels of O2, for odd-iV (a) and even-iV (b) rotational 
manifolds. The weak- field- seeking states of interest here are labeled by their Mj quantum 
numbers. 

FIG. 2. The 02( 3 S^ ) — 02( 3 S~) singlet potential for H- geometry and the contributions from 
different sets (La, Lb, L). 

FIG. 3. A Sample set of the adiabatic curves for 16 02, in this case for total angular momentum 
J = 0. In computing these curves only the values N = 1,3,5 and / = 0, 2,4, 6 are included. 

FIG. 4. Elastic scattering cross sections versus energy for different l max for 17 02- For example, 
Imax = 6 means that I = 0, 2,4, 6 partial waves were taken into account. See text for details. 

FIG. 5. Elastic partial- wave cross sections for X1 02 molecules in their magnetically trapped 
|0011, 11 > state. Odd-1 contributions do not exist for this state because of the identical- boson 
exchange symmetry. 

FIG. 6. a)Rate constants versus energy for X1 02 — 17 O2 collisions with molecules initially in 
their |0011, 11 > state. 

b) The thermally averaged elastic and loss rates from 6(a) as a function of temperature. Elastic 
collisions strongly dominate spin- changing loss collisions at low temperatures. E$ denotes the 
energy, in Kelvin units, of the height of the 17 O2 d- wave centrifugal barrier. 

FIG. 7. Elastic scattering cross sections and all inelastic scattering cross sections for the initial 
1 1122, 22 > state of 16 2 . Here channels with l max = 10, N = 1 are included. 
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